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Abstract We study the kinematical characteristics and 3D geometry of a large- 
scale coronal wave that occurred in association with the 26 April 2008 fiare-CME 
event. The wave was observed with the EUVI instruments aboard both STEREO 
spacecraft (STEREO- A and STEREO-B) with a mean speed of -240 km s-^ 
The wave is more pronounced in the eastern propagation direction, and is thus, 
better observable in STEREO-B images. From STEREO-B observations we de- 
rive two separate initiation centers for the wave, and their locations fit with the 
coronal dimming regions. Assuming a simple geometry of the wave we recon- 
struct its 3D nature from combined STEREO-A and STEREO-B observations. 
We find that the wave structure is asymmetric with an inclination towards East. 
The associated CME has a deprojected speed of ~750±50 km s~^, and shows 
a non-radial outward motion towards the East with respect to the underly- 
ing source region location. Applying the forward fitting model developed by 
Thernisien. Howard, and VourlidasI (l2006l ). we derive the CME flux rope position 



on the solar surface to be close to the dimming regions. We conclude that the 
expanding flanks of the CME most likely drive and shape the coronal wave. 

Keywords: Shock waves. Coronal Mass Ejections 



1. Introduction 

Wave-like disturbances in the solar corona were for the first time imaged by 
the Extreme-ultraviolet Imaging Telescope (EIT; iDelaboudiniere et all Il995 ) 



instrument onboard the Solar and Heliospheric Observatory (SOUO), thereafter 
called EIT waves (Moses et alx .1997 : Thompson et or, more generally. 
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EUV waves. Their generation mechanism and nature is stih an issue of debate. 
In possible scenarios it is assumed that these disturba nces are fast-mode MHD 
waves which are flare-initiated and/or CME driven (e.g. Wills-Davev and Thomps"Qnl .| 

ll999^:'Wang^.'2QQQ':'Warmuth et a/.Ll2QQll: lLong e7aZl J2QQ8l:IVeronig. Temmer. and Vrsnak[| 

[2008; Gopalswamy et al.^ 2009^), solitonsor slow- mode waves (e.g.|Wills-Davey, DeForest. and Stenflol ,| 
|2QQ7[ ). Other models suppose that these disturbances are non-wave features 
that occur due to magnetic reconfiguration associated with the erupting CME 
(e.g. iDelanned . I2QQQ| : IChen et ali I2QQ2I : lAttrill et all I2QQ9I : IPai et ali |2Q1Q[ ). 
Recently, theoretical models were developed that combine wave and non-wave 
scenarios within hybr id models (see IZhukov and Aucherd , I2QQ4I : ICohen et al . 
2QQ9l : lLiu e^"all . |2Q1Q[ ). 



The original model developed by lUchidal (| 19681 ) assumes that the 3D struc- 
ture of a shock wave in the corona sweeps over the chromosphere which would 
account for the observations of coronal as well as chromospheric wave signatures 
(Moreton waves imaged in Ha). Based on combined observations in different 
wavelengths of a coronal wave, iNarukage et all ( 2QQ2I ) found from the result- 
ing differences in propagation a 3D structure of the wave. A recent review by 



Vrsnak and Cliverl (|2QQ8f ) describes the formation of MHD waves with a 3D 



piston mechanism driven either by the CME expansion or by a flare-associated 
pressure pulse. The Solar Terrestrial Relations Observatory (STEREO) provides 
new opportunities to observe coronal waves under different vantage points. In- 
triguing results giving evidence for the 3D geometry of coronal waves are reported 
for the STEREO quadrature event from 13 February 2009. For this event the 
EUV wave could be observed from the l ateral direction as well as simultaneously 
on-disk. From these unique observations iKienreich. Temmer. and Veronid (l2QQ9h 
and Patsourakos and Vourhdasl (l2QQ9h conclude that the disturbance is a 3D 
fast-mode MHD wave which is partly driven by the flanks of the associated CME. 
They further show that the EUV wave is propagating at heights of about 100 Mm 
above the solar surface which is comparable to the coronal scale height for quiet 
Sun conditions ( Patsourakos et all . |2009[ ). In a recent study bv IVeronig et al. 



(l2010[ ) the full 3D dome of an EUV wave was identified, and interpreted to 
be freely propagating in the lateral direction after the lateral expansion of the 
CME stopped, whereas the upper part of the wave dome was permanently driven 
during the eruption. For recent reviews on the issue of EUV wave initiation 
and their nature we refer to lWills-Davev and Attrilll (|2010h . [Gallagher and Lond 
(l2010h . and lWarmuthI (l2010h . 

During the 26 April 2008 CME-flare event a coronal shock wave was imaged 
with the Extreme Ultraviolet Imager (EUVI) aboard both spacecraft of the 
STEREO mission, which were separated from each other by 49.5°. In this paper 
we present a study on the kinematics and geometry of the wave as seen simul- 
taneously from two different vantage points. These observations are combined 
with a forward fitting model of the erupting flux rope, which enables us to derive 
geometrical information about the 3D structure of the disturbance. 
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Figure 1. Sequence of running difference images from STEREO- A and -B observations in the 
EUVI 195 A passband. A dome shaped structure observed from STEREO-A is indicated with 
arrows. See also the onhne animation. SOLA: ms.tex; 13 January 2013; 16:53; p. 3 
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Figure 2. Traced wave fronts observed in STEREO-B (left) and STEREO- A (right) EUVI 
195 A (red lines) and 171 A (white lines) images. Black dashed circles indicate the fit to the 
first wave front as seen in 171 A. Black crosses give the estimated initiation center of the wave. 
The limb of STEREO-B is indicated in the image of STEREO-A. 



2. Data and Methods 

The EUV wave under study occurred on 26 April 2008 and was associated with 
a CME and a flare of GOES class B3.8/Ha importance SF from a source region 
located at N08E08 (Earth view). GOES soft X-ray (SXR) observations reveal 
that the thermal emission of the flare starts at 13:54 UT with a maximum at 
14:08 UT. The event could be observed from both STEREO satellites (STEREO- 
A and STEREO-B) which were separated by 49.5°. Complementary observations 
from SOHO/LASCO and EIT provide information from Earth view. We stress 
that the EUV wave under study was associated with a coronaQ and interplan- 
etary type II radio burst observed by Wind/W AVE^. The type II burst in the 
metric and interplanetary space were studied in iGopalswamv et all ( 2009b[ ). 

The EUVI instrument and the coronagraphs CORl and C0R2 are part of 
the Sun Earth Con nection Coronal and Heliospheric Investigatio n (SECCHI; 
Howard et a/.U2008[ ) instrument suite onboard the STEREO mission (IKaiser et al 
20081 ). launched in October 2006. The STEREO mission consists of two identical 
spacecraft, which orbit the Sun ahead (STEREO-A) and behind (STEREO-B) 
the Earth near the ecliptic plane, and slowly separate from each other by ~44° 
per year. EUVI observes the chromosphere and low corona in four differen t EUV 
emission lines out t o 1.7 Rs (with Rg the solar radius) (|Wuelser et all 12004 : 



Howard et all l2008f ). During the coronal wave event under study, the EUVI 



imaging cadence was 2.5 min in the 171 A and 10 mm m the 195 A filter. For 
the comparative analysis of the early evolution of the associated CME in the 
low corona we use data from the STEREO/SECCHI inner coronagraph CORl, 



^ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SOLAR_RADIO/BURSTS/| 
^http://cdaw.gsfc. nasa.gov/CME_list/radio/waves-type2. html I 
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a classic Lyot internally occulting refractive coronagraph with a field-of-view 
(FOV) from 1.4 to 4 7? ^: and C0R2 t he ex ternally occulted coronagraph with 
a FOV of 2 to 15 Rs ( Howard et all |2QQ8[ ). For our study we used polarized 
brightness C0R2 images. In addition, SOHO /LASCO C2 (FOV: 1.5 R,-6Rs) 
and C3 (FOV: 3.5Rs-30Rs) data are studied (|Brueckner et a/.LIl995h . 



The coronal wave is tracked by manually detecting the wave fronts, separately 
for SECCHI/EUVI-A and -B observations in both the 195 A and 171 A passband, 
respectively. From a circular fit to the fir st wave front, the in itiation center of 
the wave is derived (for more details see IVeronig et all |2QQ6[ ). Figure [2] shows 
direct images of EUVI-A and EUVI-B in the 195 A passband together with the 
wave fronts determined from EUVI 171 A and 195 A running difference images. 
The distance of the wave to the initiation center is calculated by averaging the 
measured distance of each point of the tracked wave front from the derived 
initiation center along the spherical solar surface. 

In order to infer the direction of mot ion of the CME we apply the triangulation 
method developed by Temmer. Preiss. and Veronig (|2QQ9h) on the distance-time 
measurements of the leading edge of the CME observed from LASCO, STEREO- 
A, and -B. This method uses the information of the spacecraft separation angles 
and takes the propagation direction of the CME as a free parameter, in order 
to transform SOHO/LASCO distance-time measurements to STEREO- A and 
STEREO-B view. By comparing the transformed distances with the distances 
actually observed from STEREO, the direction of motion, which gives the min- 
imum deviation, results in the best estimation of the true CME propagation 
direction. To obtain an independent result for the propagation direction of the 
CME as well as to derive the surface location of the CME flux rope we also use the 
flux ro pe fo rward fltting model developed bylThernisien, Howard, and Vourlidas 
( 2006 ) and iThernisien. Vourlidas, and Howard! (|2009[ ). This model is a ravtrace 
simulation method which enables us to compute synthetic total and polarized 
brightness images using the Thomson scattering formulae from an assumed 
electron density model. Since the graduated cylindrical shell (GCS) model is 
a reasonable simulation of a flux-rope CME it can be applied to investigate 
the appearance of a CME. Characteristic CME parameters (width, cone angle, 
latitude, longitude) are derived by fitting the density model until a best match 
is found for contemporaneous image pairs from STEREO- A and STEREO-B, 
which observe the CME from two different vantage points. 



3. Results 

Figure [T] shows a sequence of EUVI 195 A running difference images from 
STEREO- A and -B (see also the online movie). The location of the AR from 
which the wave is launched is EOS (if not stated otherwise, heliographic coor- 
dinates refer to Earth view). Hence, from the vantage point of STEREO-B the 
region is located less than 15° west from Sun center, whereas from the vantage 
point of STEREO- A it is more than 35° East. Therefore, measuring the wave 
fronts from STEREO-B is less affected by p rojection effects originating from , 
the three-dimensional nature of the wave (see Kienreich, Temmer, and Veronigl , 
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Figure 3. Wave kinematics for the western direction as derived from EUVI 171 A and 
195 A images of STEREO- A and STEREO-B (cf. Figure [2]). Solid and dashed hnes show 
the quadratic and hnear fit, respectively, to the data points. The mean velocity, resulting 
from the linear fit is given in the legend. 



2QQ9I : IPatsQurakos et all I2QQ9I ). Both STEREO- A and -B observations reveal 



that the wave intensity is higher for the eastern propagation direction. The wave 
appears to be brighter as weh as more diffuse in 195 A compared to the 171 A 
passband. In general, coronal dist urbances are found to be be t ter ob served in 
195 A than in other wavelengths ( Wills-Davev and Thompson! , 19991 ). We also 
note that a dome shaped structure is identified in STEREO- A observations which 
images the disturbance from an almost lateral direction. 

Figure [2] shows all tracked wave fronts for each spacecraft extracted from 
images in the EUVI 195 A as well as 171 A filter. From the vantage point of 
STEREO-A the AR is close to the eastern limb, therefore, the wave is best 
observed for the western direction (although less intense) and we focused in the 
construction of the wave kinematics to this direction. STEREO-B imagery enable 
us to track the wave in both the eastern and western propagation direction. 
By fitting circles to the first wave fronts separately for the eastern and western 
propagation direction observed in STEREO-B, two initiation centers on opposite 
sides of the AR are derived. This is obtained from EUVI-B 195 A as well as 171 A 
observations. 

Figure [3] shows the derived wave kinematics for the western propagation 
direction obtained from STEREO-A and -B. In the 171 A passband the wave 
can be tracked over a shorter distance than in 195 A due to the weaker signal. 
For this propagation direction distances as derived from the 195 A filter devi- 
ate from those measured in 171 A. This is found from STEREO-A as well as 
STEREO-B measurements. Especially the first front observed in STEREO-A 
195 A differs clearly from a circular shape since part of the wave is detected off 
limb. Performing a linear fit over all measurements we obtain a mean velocity of 
-245±20 km s-^ The (point-like) center of initiation from STEREO-A 171 A 
observations is derived at [-544'',193''], from STEREO-B 171 A at [318'',306''] 
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Figure 4. Wave kinematics for the eastern direction as derived from EUVI 171 A and 
195 A images observed with STEREO-B. Sohd and dashed hnes show the quadratic and 
hnear fit, respectively, to the data points. The mean velocity, v, resulting from the linear fit 
to the distance-time measurements is given in the legend. 



from Sun center. A transformation of the coordinates from A to B view gives 
[253'', 258''], i.e. there is a mismatch of ^bO" when deriving the wave initiation 
center from the different vantage points. We assume that this difference is due 
to projection effects, uncertainties in the manual tracking of the wave as well 
as due to the sini plified assumption of a point-like center of initiation (see e.g. 
Muhr et allhOld) . 

Figure m shows the derived wave kinematics from STEREO-B observations for 
the eastern propagation direction. The mean velocity of the wave is ~240±12 km s~^| 
which is comparable to the western propagation direction. The wave extracted 
from 171 A images can be tracked up to the same distance as in the 195 A 
passband. Since for the western propagation direction we observe the wave well 
in STEREO-A, we can exclude that the wave "disappears" behind the STEREO- 
B limb. Therefore we infer an asymmetry in the intensity of the wave, i.e. the 
eastern propagation direction is more intense. In the following we will show that 
the main direction of the CME propagation as well as the 3D structure of the 
wave are in accordance with the more intense eastern wave propagation direction. 



3.1. Propagation direction of the associated CME 

In Figure [5] the result of the forward fit to the STER EO-A and -B CORl image 

pair re cord ed at 14:45 UT is shown using the model by Thernisien, Howard, and Vourlidcis^ 
(|2QQ6h and iThernisien. Vourlidas. and Howard! (|2009|). The CME observations 
from 14:45 UT do not overlap in time with the EUV wave, but are only ~20 min 
after the last detected wave front. CORl observations from STEREO-A and 
-B would show the CME in earlier images. However, the more of the CME 
structure is visible the more reliable is the forward fit. As can be seen in the 
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Figure 5. Flux rope (green mesh) resulting from the raytracing- fitting method applied to 
CORl observations from STEREO-A (right) and STEREO-B (left) overplotted on composite 
EUVI 195 A -CORl difference images from ~14:45 UT. Arrows outline the shape of the 
CME to which the flux rope was fitted. See also the online animations of STEREO-A and -B 
composites from EUVI and CORl running difference images. 



lower right panel of Figure [5l the structure to the right of the top arrow possibly 
corresponds to a deflected streamer and not to the CME eruption itself (see 
also the accompanying movie). Therefore, it is possible that the CME fit, which 
includes this structure, overestimates the true width of the CME. From the fit 
we derive a CME propagation direction of E20-E25 and a width of ~60° as 
well as the position of the flux rope on the solar surface (cf. Figure [7j). Figure [6] 
shows a second method from which the CM E propagation direction is der i ved us - 
ing the triangulation method developed by Temmer. Preiss. and Veronid ( 2QQ9I ). 
From the triangulation of the STEREO-A and -B CME kinematics we obtain 
a propagation direction of E28-E35 (see Figure [6]) which is consistent with the 
result derived from the flux rope forward fit. For comparison, the AR from 
which the CME is launched has a longitudinal position of EOS. Using the derived 
information of the CME propagation direction, we obtain a deprojected CME 
speed of ~750±50 km s~^. 

Figure [71 shows a base difference image from STEREO-B EUVI 195 A and 
the derived wave fronts over the time range 13:56 UT-14:16 UT together with 
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Figure 6. Triangulation method from iTemmer, Preiss, and Veronigj (|2009l ) applied to the 
derived CME height-time data from observations of different spacecraft. SoHd hnes are the 
modeled distances for STEREO-A view (red), STEREO-B view (blue), and LASCO (gray). 
The red and blue dashed lines show the resulting de-projected distances separately derived from 
LASCO/STEREO-A and LASCO/STEREO-B comparison. The separation angle between the 
STEREO spacecraft with respect to LASCO is labeled (ZAL, ZBL). 

the flux rope position on the solar surface obtained with forward fltting appHed 
to the CORl STEREO-A and -B image pairs at 14:45 UT. The EUVI base 
difference image reveals coronal dimming regions on opposite sides of the AR. 
Their location is consistent with the two initiation centers obtained for the wave. 
The position of the footpoints of the flux rope obtained from the forward fit is 
also closely located to the dimming regions. The eastern footpoint of the flux 
rope is positioned to the east from the AR and outside the dimming region. The 
western footpoint is close to the AR and within the dimming area. 

3.2. 3D geometry of the wave 

Figure [8] (top panels) shows difference images of the wave in EUVI 171 A as ob- 
served from STEREO-A and STEREO-B at 13:59 UT, respectively. For a better 
comparison the STEREO-B image is rotated to the viewing angle of STEREO- 
A. Imaging the wave from two different vantage points reveals differences in the 
morphology of the wave as well as "foreshortening effects" . This hints at the 3D 
nature of the wave, and can be used to infer information of its 3D geometry. For 
the sake of simplicity we assume that the wave consists of a base and two lateral 
borders which are inclined with respect to the vertical. As can be seen from the 
red bars in Figure [8] the distances between the active region and the eastern and 
western part of the wave are similar observed from STEREO-B, but different 
observed from STEREO-A. We therefore may assume that STEREO-B observes 
the base of the wave nearly free from projection effects. Using the observed 
extension of the eastern and western part of the wave as seen from STEREO-A, 
and knowing the location of the active region as well as the spacecraft positions, 
we can simply use trigonometric functions to derive the inclination of the borders 
of the wave. The results are schematically drawn in the bottom panel of Figure [H 
revealing that the wave structure is asymmetric, i.e. more inclined towards East. 
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Figure 7. Traced wave fronts (13:56UT-14:16UT) from STEREO-B EUVI 195A in blue and 
17lA in white plotted on a base difference image. Red dashed circles give the fit to the first 
wave front observed in EUVI 171 A from which the initiation centers (red crosses) of the wave 
are obtained within the bipolar dimming regions. Results from the raytracing tool derived 
from CORl images at 14:45 UT (cf. Figure O are indicated in yellow: shaded circles mark the 
position and extension of the footpoints of the fiux rope, cross and dashed line give the apex 
and inclination of the fiux rope projected on the solar surface. 



4. Discussion and Conclusion 

For the EUV wave that occurred on 26 April 2008 we obtained two initiation 
centers located on either side of the AR from which the wave was launched. The 
two initiation cente rs are located within the coronal dimming regions. A similar 
result was found bv lMuhr et all ([2010: ) for the Ha Moreton wave of 28 October 
2003. For this event the initial location of the wave center was derived to be 
located close to the dimming regions well outside from t he AR core, which are 



indic ative of the footpoints of the flux rope of a CME fe.g. ISterling and HudsonI , 
19971 ). Model simulations of fast-mode MHD waves by IWand tOOi£ showed 



that the wave disturbances needed to be launched in the periphery of the AR 
to expand horizontally over the surface, whereas the initiation in the core of 
the AR with strong magnetic fields would lead to a vertical upward motion of 
the wave. These results are in accordance with the observations and the ob- 
tained initiation centers on opposite sites of the active region in weak magnetic 
field. Simultaneous imaging from t wo different vantage points enabled us to 
appl y the forward fitting model by Thernisien. Howard, and VourlidasI ( 20061 ) 



and Thernisien. Vourlidas. and Howard! ( 2009 ). From this we derived the fiux 



rope position on the solar surface and found that it is close to the coronal 
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Figure 8. Top: EUVI-A and -B 17lA difference images of the wave at 13:59 UT. Note that 
the STEREO-B image (left) is ahgned to the viewing direction of STEREO-A (right) in order 
to make the differences in the morphology of the wave as viewed from two different vantage 
points clearly visible. The lateral extension of the wave is schematically indicated with red 
bars. Bottom: Reconstructed 3D geometry (view from solar north pole) with a ~ 162° and 
(3 ~ 111°. AR indicates the location of the active region. For details of the reconstruction see 
text. 
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dimming regions and wave initiation centers, respectively. Theses findings sug- 
^e st that the CME expanding flanks are initia ting the wave. A recent study 
bv IPatsourakos. Vourhdas. and StenborgI ( 2Q1Q[ ) showed that CMEs undergo a 
strong but short Hved lateral overexpansion in their early evolution. 

The mean speed of the wave of ~240 km s~^, which is in the range of the 
characteristic speed of the ambient quiet solar corona, is derived to be similar 
for the western and eastern propagation direction as well as from observations of 
two different vantage points. The wave shows no clear deceleration which is due 
to its slow initial speed and, as a consequence of the extreme solar minimum, 
due to t he weak magnet i c field environment in which the wave is propagating 
(see e.g. iKienreich et all 12011 ). The derived speed of the wave lies within the 
velocity range of of 210-350 km s~^ as expec t ed fo r fast magnetosonic waves 
during quiet Sun conditions (see iMann et all I1999I ). The radial speed of the 
upward moving CME is much larger and of ^750 km s~ ^. In this respect it 
would be interesting to do a comparison to the model by lAttrill et all ( 2007 ) 
who interprets EUV waves as magnetic footprints of CMEs. This model would 
require a motion of the EUV wave similar to the global average motion defined 
by the lateral expansion of the CME. However, the CME is observed off-limb 
well after the end of the EUV wave and therefore it is not possible to compare 
simultaneous measurements of the wave and CME widths. 

For the western propagation direction, the measurements fr om the 195 A pass- 
band show systematic larger distances as compared to 171 A. IPatsourakos et all 
( 2OO9I ) also found differences in the wave appearance between 171 A and 195 A, 
especially in the beginning of the wave evolution. For the eastern propagation 
direction the wave is observed to be more intense. The transformation of EUVI 
images from STEREO-B view to STEREO-A view makes foreshortening effects 
clearly visible (see Figure [8]). Reconstructing the 3D structure of the wave by 
applying simple geometry, we find that the wave is inclined towards East. These 
results indicate that the 3D structure of th e wave is not symmetric. From two 
i ndepe ndent methods (forwa rd fitting model iTh ernisien. Howar d, and Vourlidas 
(|2006[ ) and triangulation Te mmer. Preiss. and Veronig (2009)) we derived that 
the propagation direction of the associated CME deviates eastward from the 
AR position. The appearance of the wave, especially its asymmetric intensity 
and 3D geometry, seems to be closely related to the non-radial evolution of the 
associated CME. The results from the present study support the wave nature 
of the phenomenon as well as the close association of the wave structure to the 
evolutionary characteristics of the expanding CME fianks. 

Since in the EUVI 171 A passband wave featu res might be confused with 
expanding loops (^Patsourakos and Vourlidasl . [2QQ9l ). we would like to emphasize 
that the double-initiation center is found also from wave fronts identified in 
195 A images. Moreover, the smooth kinematics observed in the eastern as well as 
western propagation direction gives no indication that the first fronts identified 
would not be associated to the propagating wave. In this respect we would also 
like to note that, if the first measured fronts are expanding loops, we face a 
very smooth transition between loop expansion and the build up of the wave. 
We expect that detailed studies with high spatial and temporal resolution SDO 
data will shed more light on the initiation phase of CMEs and their association 
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to coronal waves. [lIu et al. ( 2Q1Q[ ) presented the first results from SDO imagery 
for the 8 April 2010 EUV wave and proposed a hybrid hypothesis combining 
both wave and non-wave aspects to explain the observations. 
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